Astronomy & Astrophysics manuscript no. fclli's235 


February 2, 2008 


(DOI: will be inserted by hand later) 





New light on the S235A-B star forming region* 

■ Marcello Felli 1 , Fabrizio Massi 1 , Alessandro Navarrini 2,3 , Roberto Neri 3 , Riccardo Cesaroni 1 , Tim Jenness 4 

o 
o 

' 1 INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5, 50125 Firenze, Italy 

, | 2 Radio Astronomy Laboratory, University of California, Berkeley 601 Campbell Hall, Berkeley, CA 94720, USA 

3 IRAM, 300 Rue de la Piscine, Domaine Universitaire, F-38406 St. Martin d'Heres Cedex, France 

4 Joint Astronomy Centre, 660 N. A'ohoku Place, Hilo, HI 96720, USA 



Received; Accepted 



> 

m 
oo 



Abstract. The S235A-B star forming region has been extensively observed in the past from the radio to the 
near IR, but what was happening in the immediate surroundings of the water maser, placed in between the two 
nebulosities, was still unclear because of insufficient resolution especially in the spectral range from the Far IR 
to the mm, even though there were sound indications that new young stellar objects (YSOs) are being formed 
there. We present here new high resolution maps at mm wavelengths in different molecules (HCO + , C 34 S, H2CS, 
£f~) • SO2 and CH3CN), as well as in the 1.2 and 3.3 mm continuum obtained with the Plateau de Bure interferometer, 

and JCMT observations at 450 /im and 850 fim that unambiguously reveal the presence of new YSOs placed in 
between the two Hn regions S235A and S235B and associated with the water maser. A molecular core and an 
unresolved source in the mm and in the sub-mm are centred on the maser, with indication of mass infall onto 
the core. Two molecular bipolar outflows and a jet originate from the same position. Weak evidence is found 
for a molecular rotating disk perpendicular to the direction of the main bipolar outflow. The derived parameters 
indicate that one of the YSOs is an intermediate luminosity object (L ~ 10 3 Lq) in a very early evolutionary 
phase, embedded in a molecular core of ~ 100 M©, with a temperature of 30 K. The main source of energy for 
the YSO could come from gravitational infall, thus making of this YSO a rare example of intermediate luminosity 
protostar representing a link between the earliest evolutionary phases of massive stars and low mass protostars 
of class 0-1. 
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1. Introduction et al. 1990), a well known landmark of recent star for- 

_„__ . , Mrio _ , . „ mation (see e.g. Tofani et al. |[lj?951). With the resolution 

S235A and S235B are two close-by (separation ~r) small ,, , , ,,,,,, ., , , , 

„,/,:, available at that time it was not clear whether the maser 

patches of nebulosity (—20 and —10 m diameter, respec- . , , .,, „„,,., c not . n 1 1 i , 

. , , , ,1 „» , was associated with S235A, S235B or simply happened to 

tively), so called because they are located about 10 south , ,, , r -/tip t? n- n iwm 

, , , , TT . , mrwm De 111 the same star forming region (Israel & Fclli 1978). 

ol the more extended Hn region S235 (bharpless 1959. c . , , , , , , , , , c 

, ° , ,i i- , Single dish observations have also revealed the presence oi 

Felli & Lhurchwell 1972), but most probably not directly ,, , / AT , „ , n A6^ i tt i • i x 

. /T », „ - . ' ; a methanol maser (Makano & Yoshida 1986; Haschick et 

related to it (Israel & Felli 1978). About 3.5 arcmm south , nm j c a-r\ ru • t, i rnrrrai 

, V 1 ' al. 1990) and of an biO maser (Harju et al. 1998). 

of S235A-B there is another patch of nebulosity called 

S235C and identified with a small Hn region (Israel & To ascertain the true relationship between S235A-B 
Felli [UT75|l . and the water maser and establish which is the source of 
Early interferometric radio continuum observations energy supply for the water maser emission, high resolu- 
have shown that S235A is a young Hn region, even though tion radio observations were performed with the VLA in 
far beyond the earliest ultracompact (UC) evolutionary the A configuratio n in t he water line and in the contin- 
phase, while no radio continuum emission comes from uum (Tofani et al. 11995(1. which led to the result that the 
S235B (Israel & Felli fTWHjl . The interest in these two neb- H 2 maser is well outside the boundary of the S235A Hn 
ulosities was increased by the fact that an H 2 maser is region, away from S235B and that no point-like radio con- 
present in the same area (Henkel et al. Hrjgpl Comoretto tinuum source (down to a flux density of -0.5 mJy and 

with a resolution of C/3) is associated with it, i.e. no trace 

Send offprint requests to: M. Felli, mfelli@arcetri.astro.it q£ ft uc Hn region near the maser wag found 

* Based on observations carried out with the IRAM Plateau 

de Bure Interferometer. IRAM is supported by INSU/CNRS At that time it was becoming clear that water masers 

(France), MPG (Germany) and IGN (Spain) without associated compact radio continuum emission 
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could trace one of the earliest phases in the formation 
of a young stellar object (YSO), in which the UC Hn re- 
gion ionized by the early type star (if at all present) is 
so dense and optically thick to escape detection at radio 
wavelengths (Codella et al. H^MlHM7|l . 

The presence of a molecular cloud in the S235A-B re- 
gion was known from earlier low resolution observations 
in several molecular lines (Evans & Blair ITMT1 Lafon et 
al. 113531 Nakano & Yoshida H55S|l . To test the hypothesis 
that the water maser indicated the location of one or more 
new YSOs independent from S235A-B, and get more in- 
formation on the environment in which the water maser 
was located, new maps in molecular tracers of high den- 
sity gas (Cesaroni et al. I1999b() and near-IR observations 
of the region (Felli et al. H997fl were made. 

The results confirmed the validity of the hypothesis: 
the water maser is located at the centre of a high density 
small diameter molecular core, which peaks just between 
the S235A and S235B nebulosities. A cluster of near-IR 
sources is also located in the same area. The source with 
the largest PR excess (called Ml) is placed close to the 
water maser. The picture before the present observations 
is summarized in Fig. which shows the overlay of the 
HCN map obtained with a resolution of 27" (Cesaroni 
et al. I1999b|l with a K-band (2.2 /j,m) map of the same 
area (Felli et al. I1997f) . All the indications suggest that 
the water maser, Ml (both indicated in Fig. ^) and the 
HCN clump are all part of a much younger episode of star 
formation with respect to the S235A-B nebulosities, not 
necessarily directly related to them. 

For completeness it should be noted that the VLA ob- 
servations of the water maser (Tofani et al. 11995(1 refer to 
the component at ~ —60 km s _1 , which happened to be 
flaring at the time of the observations (1992) and was the 
only feature above the noise in the observed bandwidth 
(from -78 to -40 km s _1 ). During the continuing patrol 
of the water maser emission with the Medicina radiotele- 
scope, started in 1988, this velocity component has very 
seldom emerged again from the noise level, while two other 
velocity components between ~ —2.5 and 2.5 km s^ 1 are 
frequently observed (Felli et al., in preparation). A partial 
report of the results of the Medicina patrol is given in Felli 
et al. ifTMTjl . 

In the present work we investigate with better resolu- 
tion the region around the water maser to study in more 
detail the YSOs that presumably are harboured in the 
centre of the molecular cloud. We present arcsec resolution 
maps at mm wavelengths in different molecules that trace 
high density gas (HCO+, C 34 S, H 2 CS, S0 2 and CH 3 CN), 
as well as in the 1.2 and 3.3 mm continuum sensitive to 
thermal dust emission. 

To obtain the bolomctric luminosities of the YSOs and 
their spectral energy distribution (SED), high resolution 
far IR observations are needed, since at these wavelengths 
YSO spectra are expected to have their maximum. We 
inspected the IRAS HIRES maps (Aumann et al. QW01 
of the region around the water maser (which contains 
IRAS 05375+3540) in the four IRAS bands, but the res- 




c<(j2G0O) 

Fig. 1. Overlay of the Pico Veleta HCN map (contours, 
resolution 27") from Cesaroni et al. (1999b) with the 
TIRGO K-band image of the S235A-B region (grey scale) 
from Felli et al. (|1997[> . The VLA position of the wa- 
ter maser from Tofani et al. ( 1995) is indicated with a 
cross (much larger than the positional uncertainty which 
is 0.1"). Ml, the K-band source with the largest infrared 
excess, is also indicated. 

olution is inadequate for separating the contribution of 
S235A and S235B from that of the YSO. 

We also present here archive JCMT observations at 
450 /xm and 850 [im that unambiguously demonstrate that 
the sub-mm emission peaks on the YSO associated with 
the water maser and separate its emission from that of 
S235A-B. 

For consistency with previous work a distance of 1.8 
kpc will be assumed (Nakano & Yoshida 1986). 

2. Observations and results 

2.1. Plateau de Bure observations 

The interferometric observations were carried out on 
November 22, 1997 and March 25, 1998 in the CD set 
of configurations of the IRAM five element array at 
Plateau de Bure, France (Guilloteau et al. I1992J) . The 
five 15-m antennas were equipped with 82-116 GHz and 
210-245 GHz SIS receivers operating simultaneously with 
double-sideband (DSB) noise temperatures of 35 K at 3.3 
mm and 50 K at 1.2 mm. The receivers were tuned DSB 
at both wavelengths. The facility correlator was centred 
at 89.014 GHz in the lower side-band at 3.3 mm and at 
241.102 GHz in the lower side-band at 1.2 mm. We show 
the frequency set-up of the correlator and list the detected 
spectral lines in Table ^ The six units in the correlator 
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were placed in such a way that a frequency range free of 
lines could be used to measure the continuum flux. The ef- 
fective spectral resolution is about twice the channel spac- 
ings given in Table ^ Source visibilities were phase cali- 
brated by means of interspersed observations (every 20 
minutes) of nearby point sources. 

The bandpass, amplitude and phase calibration was 
carried out in the standard antenna based manner. The 
flux of the primary calibrator was bootstrapped from 
IRAM monitoring measurements and used to derive the 
absolute flux density scale. Table[2]gives a list of the main 
parameters for our Plateau de Bure interferometer obser- 
vations. The phase centre position is that of the water 
maser (Tofani el al. 0.995). 

The calibration and data reduction were made using 
the GAG software developed at IRAM and Observatoire 
de Grenoble. Continuum subtraction was performed in the 
UV plane by using the integral over the line-free channels 
of the 160 MHz units. Finally, channel maps were pro- 
duced for all the lines. The conversion factor from flux 
to brightness temperature in the synthesized beam is ~ 
5.9 K (Jy/beam) _1 at 3.3 mm and ~ 5.8 K (Jy/beam) -1 
at 1.2 mm. 

2.1.1. Continuum 

The 3.3 mm and 1.2 mm continuum maps (uncorrected 
for primary beam attenuation) are shown in Fig. [2J 

In the 3.3 mm map, S235A is clearly detected as an 
extended source at the edge of the primary beam. The 
total flux density (corrected for primary beam attenua- 
tion) is 54 mjy This value is lower than the optically thin 
extrapolation of the free-free emission observed at lower 
frequencies (^200 mjy using the flux densities of Israel & 
Felli 11978")) . but confirms that the radio continuum emis- 
sion comes from the ionized gas of the S235A Hn region, 
with negligible contribution from dust emission. The rea- 
son for the lower observed flux density must probably be 
searched for in the higher resolution of the present ob- 
servations (5'.'6 compared with 12") which misses lower 
surface brightness extended emission. The overlay with 
the K-band image (Felli et al. 11997)1 . not shown here, in- 
dicates that the S235A mm emission is centred on the 
western edge of a K-band nebulosity, where a sharp drop 
in surface brightness suggests the presence of a foreground 
dust front. 

S235B is not detected at 3.3 mm with an upper limit 
to a point source at the same position of 1 mjy and also 
undetected at 1.2 mm, but in this case it is also well 
outside the primary beam. The nondetection is consis- 
tent with previous upper limits at lower frequencies (Israel 
& Felli I1978|l and confirms that no radio emission either 
from ionized gas or dust is coming from S235B. This be- 
haviour of S235B at radio wavelengths is surprising, since 
the source has been detected in Ha (Glushkov et al.[1975) 
and Br7 (Krassner et al. 1982). Possible explanations for 



this effect in terms of an ionization-bounded expanding 
envelope have already been discussed by Felli et al. 0997). 

S235C is totally outside the primary beam both at 3.3 
and 1.2 mm. 

The new important result of the present observations, 
in connection with the hypothesis of YSOs associated with 
the water maser, is the detection at 3.3 and 1.2 mm of 
an almost unresolved core slightly to the north (2") of 
the water maser. At the highest resolution of the 1.2 mm 
observations the core diameter is 3", only slightly larger 
than the synthesized beam. There is also a tail protruding 
approximately to the south of the core, more evident in the 
3.3 mm map, but also partly visible in the 1.2 mm map, 
extending for ~10" in the longest direction and almost 
unresolved in the other. 

Because of the core-tail morphology we have tried to 
separate the flux density of the core from that of the tail. 
The flux densities of the core (corrected for primary beam 
attenuation) are 20 and 245 mjy at 3.3 and 1.2 mm, re- 
spectively. The spectral index between 3.3 and 1.2 mm 
is 2.5. Since we cannot estimate the possible contribution 
to the 3.3 mm continuum by free-free emission, this is a 
lower limit to the dust spectral index, but it clearly indi- 
cates that most of the emission comes from the dust that 
surrounds the YSOs. The extrapolated flux density of the 
core dust emission at 8.4 GHz is 0.06 mjy, well below the 
sensitivity of the VLA observations of Tofani et al. (1995). 
Thus the contribution of dust emission at 8.4 GHz is likely 
to be negligible with respect to free-free emission, so that 
the upper limit of 0.5 mjy (Tofani et al. H995|l can still be 
used to constrain the emission of a UC Hn region at the 
position of the water maser. 

For the tail, the flux densities are 7.0 and 13 mjy at 3.3 
and 1.2 mm, respectively, with some uncertainty in the 1.2 
mm flux density because of the weakness of the emission. 
The spectral index of the tail between 3.3 and 1.2 mm is 
0.6, definitely different from that of the core and identi- 
cal to that expected either from an ionized stellar wind 
(Panagia & Felli 11975)1 or a thermal jet (Reynolds , 
and suggests that the emission comes from an ionized jet, 
presumably departing from the core, with very small or 
null contribution from dust emission. As a counter-check, 
the predicted flux density of the jet at 8.4 GHz using a 
spectral index 0.6 is of the order of 1.7 mjy, so the fact 
that it was not detected in previous observations at lower 
frequencies is consistent with their sensitivities and reso- 
lution. Deeper radio continuum observations at cm wave- 
lengths are essential to better define the parameters of the 
thermal jet. 

2.1.2. HCO+ 

The HCO + (1-0) map of the emission averaged over the 
entire velocity range (from -49 to 11.5 km s _1 ) is shown 
in Fig. The line emission peaks at the position of the 
water maser, making tighter the connection between water 
maser and molecular core than what was found at lower 
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Table 1. Frequency set-up for the Plateau de Bure interferometer. 



Line 


Frequency 


Band 


Channel 


spacing 11 




(MHz) ' 


(MHz) 


(MHz) 


(km s" 1 ) 


3.3 mm continuum 


90494 


160 


2.5 


8.29 


HCO+(1-0) 


89189 


20 


0.078125 


0.26 


CH 3 CN(5-4) 


91987 


40 


0.15625 


0.51 


1.2 mm continuum 6 


242631 


160 


2.5 


3.09 


C 34 S(5-4) 


241016 


40 


0.15625 


0.19 


H 2 CS 7(l-6)-6(l,5) 


244048 


160 


2.5 


3.07 


S0 2 14(0,14)-13(1,13) 


244254 


40 


0.15625 


0.192 



a the effective spectral resolution is about twice the nominal channel spacing 

6 the 1.2 mm continuum was obtained by averaging two 160 MHz bandwidths in the lower and upper side bands 



Table 2. Instrumental parameters for the IRAM Plateau de Bure interferometer observations 



Parameter 



Value 

= 05 h 40 m 53?420 



Centre of phase (water maser) 

Observing mode 
Number of antennas 
Baseline range 
Band centre 

Number of sections in the correlator 
Primary HPBW 
Synthesised HPBW (PA.) 

Primary flux density and 
bandpass calibrators 

Phase calibrators 



a (2000) 
(5(2000) = 35°41'48'.'80 
continuum+line; double side-band at 3.3 mm and at 1.2 mm 

5 

15-150 m 

89.014 GHz (LSB); 241.102 GHz (LSB) 
6 

55"; 20" 
3.3 mm: 5" 4 x 4'.'2 (61°) 
1.2 mm: 2'.'2 x l'.'8 (45°) 
0415+379 (5.5-6.7 Jy at 3.3 mm; 3.4-4.1 Jy at 1.2 mm) 
3C273 (25.3-29.8 Jy at 3.3 mm; 13.2-26.2 Jy at 1.2 mm) 

3C345 (4.1 Jy at 3.3 mm; 1.1 Jy at 1.2 mm) 
2013+370 (1.8-2.2 Jy at 3.3 mm; 0.4-1.3 Jy at 1.2 mm) 
2005+403 (0.65-0.75 Jy at 3.3 mm; 0.08-0.24 Jy at 1.2 mm) 



resolution in other molecular lines (Cesaroni et al. 1199 9b). 
The molecular core has a FWHM diameter of ~ 11" and 
is surrounded by a more extended lower intensity region 
with a FWHM size of - 30". 

Comparison of the HCO + and the 3.3 mm continuum 
maps does not show any obvious evidence of connection 
between the molecular core and S235A (e.g. no sharp rim 
either in the molecular cloud or in the Hn region at the 
location where the two are closest), thus making it very 
unlikely that the ionized gas of S235A and the HCO + 
molecular cloud might be interacting. 

Inspection of the line profiles at two different positions 
near the water maser (see Fig. reveals three distinct 
features: 1) extended wings suggesting the presence of an 
outflow; 2) a shift in the velocity of the peak with posi- 
tion, which suggests the presence of two distinct molecular 
components at different velocities; and 3) a dip of emission 
on the red side of the profile at -15.5 km s _1 , suggestive 
of red-shifted self-absorption and hence of infall inside the 
core. 

To study the morphology of the outflow we have in- 
spected carefully the channel maps and produced a map 
averaging all channels up to the edge of the band with 
w lsr > -15.8 km s _1 for the red lobe and with ^lsr < ~~ 



20.8 km s _1 for the blue lobe (see full lines in Fig.^J. The 
result is shown in Fig. El Two bipolar outflows are evident: 
1) a more extended one oriented in the NE-SW direction 
and 2) a smaller one oriented in the NNW-SSE direction. 
Within the uncertainty on the separation of the two out- 
flows, both seem to originate from the position of the water 
maser. The flux contained in the NE-SW outflow lobes is 
~12% of that over the entire velocity range. The presence 
of a larger scale bipolar outflow in the S235A-B region 
had also been reported by Nakano & Yoshida fl^86 ) using 
the 12 CO line and by Cesaroni et al. I|1999b|l and Felli et 
al. I|1997[) with 13 CO observations. However there is very 
little correspondence between those maps and the map 
of Fig. most probably the different resolutions and the 
fact that the CO line is much more affected by the large 
scale, lower density molecular gas make it impossible to 
compare the results of Cesaroni et al. ( 19 99b|) and Felli et 
al. I|1997|) to the map in Fig. El 

To investigate the origin of the emission at different 
velocities (see Fig.0J, we produced two maps in two small 
velocity intervals centred around the velocities marked 
with dotted lines in Fig. 01 (approximately -17 and - 
19 km s _1 ). The two maps are shown in Fig. |S^,-b and 
cover the velocity intervals from -18.16 to -15.8 km s _1 
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5 h 40 m 54 s 5 h 40 m 52 s 

a(J2000) 

Fig. 2. The 3.3 mm (a) and 1.2 mm (b) maps. The water maser position is indicated by a cross. The dotted circle 
defines the primary beam HPW. The box in the lower left shows the synthesized HPBW. The position of S235A 
(detected) and S235B (undetected) are marked in the 3.3 mm map. Contour levels range from 1 (3 a) to 9 by 1 
mjy/beam at 3.3 mm and from 3 (3 a) to 57 by 9 mJy/beam at 1.2 mm. The dashed square box in the 3.3 mm map 
delimits the same area displayed in the 1.2 mm map. 



and from -20.8 to -18.16 km s _1 ; each of them presents a 
peak, but at a different position, suggesting the presence 
of two distinct cores which will be referred to as C17 and 
C19. The separation between C17 and C19 is ~7" in the 
NS direction, with C17 (the most intense) slightly to the 
north of the water maser position and C19 (the weaker 
one) to the south of it. 

Finally, to establish if the dip on the red side of the 
profile is due to infalling material, we have produced two 
maps from channels symmetrically placed in velocity with 
respect to the peak HCO + velocity (-17.2 km s" 1 ), i.e. at 
-15.5 km s _1 (the velocity of the dip) and -18.9 km s _1 
(the symmetric velocity). The two maps are shown in 
Fig. in grey scale, overlaid on the CH3CN map (con- 
tours), which traces the densest part of the molecular core. 
The dip clearly coincides with the peak of the CH3CN 
emission, while at the symmetric velocity there is no de- 
crease of HCO + emission at the position of the CH3CN 
peak (apart from some complication due to the presence 
of the HCO + C19 peak discussed above). This indicates 
that the dip is a signature of infall in the core. 

The velocity of the infalling gas can be estimated us- 
ing the model of Myers et al. I|1996|l and applying their 
Eq. (9), which for the case of Td (the brightness tempera- 
ture of the dip) much smaller than Tbd (the height of the 
blue peak above the dip) and Xrx> (the height of the red 
peak above the dip) simplifies to: 



where a is the velocity dispersion of the line and (u re d _ 
fbiue) is the difference between the velocities of the red 
and blue peaks. From Fig. H T BD =2.4 K, T RD =0.12 K, 
w rc d-Vbiuc=4 km s _1 and from Fig.^]]<7=2.6 km s -1 . The 
derived V; n is ~ 5 km s _1 . 

Using Vi n and making plausible assumptions on the 
radius (i?i n ) to which V- m refers and the corresponding 
density (p — 3M/4n Rf n ), one can estimate the mass infall 
rate: 

M in = ATrRf nP V in = 3^^ (2) 

-Kin 

For M and i?j n we have used the values derived from the 
mm continuum observations, given in Table0J The derived 
value is M in ~ 10~ 2 M y _1 . 

2.1.3. CH 3 CN 

The CH 3 CN(5-4) K=l map averaged from -20.63 km s^ 1 
to -15.03 km s _1 is shown in Fig. |3t>- The diameter of 
the emitting region is ^7" after beam deconvolution. The 
two maps in the velocity intervals from -17.58 to -15.54 
km s _1 and from -19.10 to -17.58 km s _1 are shown in 
Fig. EJi-d. The CH3CN line profile obtained by averaging 
over the synthesized beam is shown in Fig. |S1 

The two bipolar outflows seen in the HCO + line are 
not detected in the CH3CN line since it is not possible 
to use the stronger K—0 and K—l transitions because of 
blending of the lines, while the higher order transitions, 



6 



M. Felli et al.: New light on the S235A-B star forming region 



T- 



t 1 1 1 1 1 1 1 r 

c - CH,CH(5-*) K— :1 V -18 —IB km^ 



35**2^8- - 



35**1'4Q" — 



3SMV26" 




C 3 7 




\ — I — | — i — I — i — | — HHH — 1 — I — | — f — I — I — | — h 



b - HC£l*{l-fl)-V-50.S l&— .13 kj-n/i 



d fciCHjCNfS-*}- k»t v ^20.5. \e> -IB km/* 



35M1 "23" 




Fig. 6. Maps of the line emission averaged in different velocity ranges. The water maser position is indicated by a 
cross. The scale is identical to that of the 3.3 mm map in Fig. [21 The dotted circle defines the primary beam HPW. 
Left panel: HCO + emission: a. C17; b. C19. Contour levels range from 0.15 to 1.8 by 0.15 Jy/beam (C17) and 0.06 
to 1.26 by 0.12 Jy/beam (C19). Right panel: CH 3 CN emission: c. C17; d. C19. Contour levels range from 0.01 to 0.05 
by 0.01 Jy/beam. 



where blending is less severe, are too weak. Also, CH3CN 
is a tracer of high density gas and not of bipolar outflows. 

As in the HCO + line the CH3CN emission shows two 
distinct velocity peaks (once allowance is made for the 
lower intensity of the line, see Fig.Et-d). Correspondingly, 
the two components C17 and C19 can be found in the 
same velocity intervals and at the same position as those 
observed in the HCO + line. 

2.1.4. C 34 S, H 2 CS and S0 2 

The maps in the 1.2 mm lines of C 34 S, H 2 CS and S0 2 
(see Table^for the frequency set-ups) are shown in Fig. El 
not corrected for primary beam attenuation. All of these 
confirm that the molecular emission peaks at the position 



of the water maser. The diameter of the core in the C S 
line is ~3'/5. 

Most of the H 2 CS emission appears to be centred at 
~ —18 km s _1 and clearly traces C17. The H 2 CS map 
shows also a tail (at a 3cr level) in the same direction and 
position as the jet observed in the continuum at 3.3 and 
1.2 mm, with emission concentrated around -19 km s _1 
and partly overlapping with C19. 

S0 2 is very faint, so no velocity structure can be de- 
rived from our data. It is concentrated towards C17, but 
it can be noted that its emission seems to be stretched in 
the NE-SW direction, roughly along the NE-SW outflow. 

The comparison between the C 34 S line profile observed 
at Pico Veleta and that at Plateau de Bure is shown in 
Fig. EH Most of the C 34 S flux comes from structures re- 
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Fig. 3. a. The HCO + map averaged over the entire ve- 
locity range. The water maser position is indicated by a 
cross. The scale is identical to that of the 3.3 mm map in 
Fig. [3 The dotted circle defines the primary beam HPW. 
The box in the lower left shows the synthesized HPBW. 
Contour levels range from 10 to 150 by 10 mJy/beam. b. 
The CH3CN(5-4) map averaged over the velocity range of 
the K=l transition. Contour levels range from 6 to 30 by 
6 mJy/beam. 



solved out in the interferometric observations. Hence, the 
molecular core detected in our high resolution observa- 
tions represents just a small (high density) part of the 
molecular cloud present in the S235A-B region. 

2.2. James Clerk Maxwell Telescope - SCUBA 

The submillimetre observations were carried out with the 
Submillimetre Common-User Bolometer Array (SCUBA; 
Holland et al. HMjU on the James Clerk Maxwell Telescope 
on Mauna Kea, Hawaii on October 24, 2000 and were 
obtained from the public archive at the CADC. SCUBA 
consists of a short wavelength array containing 91 pix- 



2 - 



1 - 



1 - 



in 



n — i — i — i — | — i — i — i — i — rr~r — i — i — i — | — i — i — i — r 
HCO + (1-0) (0.8,2) ■ 



-40 




-30 -20 

Velocity (km/s) 



Fig. 4. The HCO + line profiles, averaged over an area 
roughly equal to that of the synthesized beam, in two 
positions close to the peak. The offsets (in arcsec) with 
respect to the water maser position are given in the top 
right corner of each panel. The two dotted lines indicate 
the peak velocity at the two positions. The two full lines 
indicate the limiting velocities used to define the blue and 
red outflow lobes. 



els and a long wavelength array containing 37 pixels ar- 
ranged in a close-packed hexagon. Both arrays have ap- 
proximately the same field-of-view (2.3 arcmin) and can 
be used simultaneously. These observations used the 450- 
and 850-micron filters and were taken using the "Emerson 
II" scanning technique (Emerson^)95; Jenness et al.[2000) 
where the telescope is scanned in Nasmyth coordinates 
(in order to generate fully-sampled images) at 24"/sec 
whilst the secondary is chopped in fixed directions on the 
sky. Six chop configurations were used: 30", 44", and 68" 
in right ascension and similarly in declination. A single 
map was taken for each chop configuration, covering a 
region of 3 by 10 arcmin and took approximately 35 min- 
utes. CRL 618 was used for flux calibration (Jenness et 
al. I2P02 ) with minor corrections applied using pointing 
observations on 0552+398 straddling the science observa- 
tions. The zenith atmospheric opacity was obtained us- 
ing a fit to the CSO 225 GHz tipping radiometer data 
as described in Archibald et al. I|2002|) and was approxi- 
mately 0.07 (225 GHz). The data were reduced using the 
standard SCUBA User Reduction Facility (SURF; Jenness 
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Fig. 11. The SCUBA maps at 850 /im and 450 /im. The water maser position is indicated by a cross. The box in the 
lower right shows the HPBW. The dotted circle in the 850 /im map defines the HPW of the primary beam at 3.3 mm. 
The values of the contours are indicated in the grey wedges on top of each map. The source at the bottom of the two 
maps is a real detection and is related to the Hn region S235C (Israel & Felli H978). 



& Lightfoot lTMgi software. The half power beam widths 
(HPBWs) are 14" at 850 and 7.5" at 450 /mi. 

The SCUBA maps at 850 /im and 450 /im are shown 
in Fig. At both wavelengths the peak coincides with 
the water maser. 

To show the relationship between the mm and sub-mm 
peaks with the mid-IR emission we have overlaid in Fig. 1121 
the 450 /im and 3.3 mm maps with the 21 /im map of 
the Midcourse Space Experiment (MSX, resolution 20"). 
It is clear that the 21 /im emission comes predominantly 
from the S235A and S235B nebulosities, with very little 
contribution (if any) from the YSO associated with the 
water maser, the mm core and the 450 /im core. 

The two SCUBA maps show a compact source at the 
centre, surrounded by a more extended halo. This is more 
evident at 850 /im. To separate the two components, a ra- 
dial profile of the emission was obtained at the two wave- 
lengths by measuring the flux density contained in circular 
rings with increasing radii departing from the peak posi- 



tion. The profile at 850 /im is shown in Fig. ^| together 
with its decomposition into two gaussians of half power 
width (HPW) equal to 12" and 56". The corresponding 
flux densities for the compact source and the entire source 
are 5 Jy and 24 Jy at 850 /im and 11 Jy and 84 Jy at 450 
/im. 

The Spectral Energy Distribution (SED) from 3.3 mm 
to 8 /im is shown in Fig. 1141 In view of the different reso- 
lutions of the MSX and IRAS data with respect to the 
present observations and the results of the overlays of 
Fig. E21 the IRAS and MSX data refer essentially to the 
S235A-B nebulosities and will not be considered further. 
The fits to our data will be explained in Sect. 13. ll 

At the southern edge of the maps in Fig.^J well out- 
side the primary beam of the interferometric observations, 
another sub-mm source is present. This source coincides 
with S235C. The radio continuum flux densities of S235C 
at 21 and 6 cm are 70 and 60 mJy (Israel & Felli 115780 . 
and indicate that S235C is an Hn region, presumably 
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Fig. 5. The HCO + outflow. The water maser position 
is indicated by a cross. The scale is identical to that of 
the 3.3 mm map in Fig. [21 The dotted circle defines the 
primary beam HPW. The box in the lower left shows the 
synthesized HPBW. Contour levels range from 10 to 100 
by 10 mJy/beam (red lobe, dashed) and 6 to 60 by 6 
mJy/beam (blue lobe, full line). 



ionized by a B0. 5 star. In the same area there are also 
IRAS 05375+3536 and an MSX source. The flux densities 
at 850 /im and 450 /jm are 1.0 ± 0.1 Jy and 3.5 ± 0.5 Jy, 
respectively. This new result clearly indicates the presence 
of dust associated with S235C. In Fig.^Dwe show an over- 
lay of the 850 /im map with the 21 /jm image derived from 
MSX. The position and size of the 6 cm radio continuum 
source are marked with an ellipse. The not perfect match 
between the 850 /Ltm emission and the 21 /im and radio 
continuum emission shows that two close-by YSOs in dif- 
ferent evolutionary stages may be present in this area: 
the S235C Hn region and the 850 /im source, represent- 
ing most probably the emission from dust associated with 
a YSO in an earlier evolutionary phase. The situation is 
very similar to that occurring in the S235A-B complex. 
However, given the lack of interferometric information on 
S235C, no further discussion of this star forming region is 
possible. 

3. Discussion 

3.1. The molecular core 

One of the main results of the present work is the con- 
firmation that a high density molecular and dusty core, 
part of the larger scale molecular complex detected in the 
CO lines (see e.g the 1,5 CO observations of Cesaroni et 
al. I1999b| . is located close to the water maser. This mor- 
phology is typical of the earliest stages in the evolution of 
a YSO and proves that a new episode of star formation is 
occurring near the water maser, while other star forma- 
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Fig. 7. Gray scale maps of the HCO + emission at -15.5 
km s _1 (top) and -18.9 km s~ x (bottom). The contours 
reproduce the CH3CN of Fig. OH The water maser position 
is indicated by a cross. Gray levels are from 0.029 to 0.177 
in steps of 0.029 Jy/beam in the top panel and from 0.107 
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Fig. 8. The CH3CN line profile of the emission at the 
peak position averaged over the synthesized beam. The 
thin line is a gaussian fit to the observed lines. 
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Fig. 9. The C 34 S (a), H 2 CS (b) and S0 2 (c) maps of 
the velocity averaged emission. The water maser position 
is indicated by a cross. The scale is identical to that of 
the 1.2 mm map in Fig. [21 The dotted circle defines the 
primary beam HPW. The box in the lower left shows the 
synthesized HPBW. Contour levels range from 6 to 36 by 
6 mJy/beam for C 34 S, from 15 to 175 by 30 mJy/beam 
for H 2 CS and from 25 to 75 by 25 mJy/beam. for Sec- 



tion events like those associated with S235A and S235B 
occurred in the past and in S235A have had sufficient time 
to create a well developed Hn region. 

The temperature of the molecular core can be de- 
rived from the CH 3 CN(5-4), if=0,l,2,3 lines measured to- 
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Fig. 10. Comparison between the C 34 S profile observed at 
Pico Veleta with a resolution of 10" (thin line) and that 
observed at Plateau de Bure over a similar area with a 
resolution of ^2" (thick line) 



wards the peak position by means of the rotation diagram 
method (Boltzman plot), which assumes that all energy 
levels of CH3CN are populated according to local thermo- 
dynamic equilibrium (LTE) with a single excitation tem- 
perature. The results of the fit are shown in Fig. 1161 and 
give a temperature of ~30 K. This value is smaller than 
the 43 K peak main beam brightness temperature of the 
13 CO(2 - 1) line (Felli et al. 11997)) and suggests that the 
optically thick CO emission of the larger scale molecular 
cloud might be heated from outside by the radiation of 
the star exciting the S235A Hn region, while the molecu- 
lar core around the water maser is shielded from the Hn 
region radiation and excited only by the YSOs placed in- 
side. 



As noted in Sect. 12.1.41 the molecular core is also 
traced by H 2 CS and S0 2 , which are predominantly found 
towards hot cores (Hatchell et al. |1998 ) or the envelopes 
of massive stars (van der Tak et al. I20U3[) . In particular, 
van der Tak et al. |2003j ) found that H 2 CS arises from 
regions with excitation temperatures of ~ 50 K, which is 
not too different from our estimated temperature, whereas 
S0 2 shows a higher excitation temperature (~ 100 K). 

A large fraction of S0 2 is usually found in high- 
velocity gas (van der Tak et al. I2003f) : this aspect might 
be consistent with the NE-SW elongated shape of the 
S0 2 emission along the direction of the NE-SW outflow. 
Indeed, sulphur-bearing molecule production is believed 
to originate either from grain evaporation in the dense 
hot environment around protostars, or from shocked gas 
(Hatchell et al. H9l?gl Bachiller et aLEDOU van der Tak et 
al. I2()()3)l . From our observations with only one transition 
per molecule, we cannot state whether they mostly arise 
from warm gas or from shocked gas around the protostar, 
although part of the emission of H 2 CS is clearly related 
with the outflows (hence, with shocked gas). 



M. Felli et al.: New light on the S235A-B star forming region 



11 




5 h 40 m 55 s 5 h 40 m 50 s 

a(J2000) 

Fig. 12. Top panel: overlay of the 21 /im maps of MSX 
(contours) with the 450 /mi emission (grey scale). Bottom 
panel: overlay of the 21 /Ltm maps of MSX (contours) with 
the 3.3 mm emission (grey scale). Contours range from 
10~ 5 to 1.9 1CT 4 in steps of 2 1(T 5 Watt mT 2 sterad" 1 . 



To derive the luminosity of the YSOs and the mass of 
the molecular cloud we use the mm and sub-mm flux den- 
sities. FigureH~4lshows grey-body fits to our data, with the 
assumption of constant density and constant temperature 
in the molecular cloud. The best fit is obtained using a 
dust opacity K v oc i/ 5 and the temperature derived from 
the CH 3 CN lines. Larger dust opacities (k u oc v x ~ 2 ) pro- 
duce worse fits to the mm and sub-mm data. 

The lower dotted curve in Fig. El shows a fit to the 
mm data alone. In this fit no distinction was made be- 
tween core and jet and total flux densities were used. It 
should also be noted that the 1.2 mm observations may 
resolve out, and consequently miss, some of the extended 
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Fig. 13. The profile of the emission at 850 (j,m integrated 
over circular rings of increasing radii. The two gaussians 
(convolved with the beam) whose sum best fits the ob- 
served profile are represented by the dotted lines. The 
vertical dashed line represents the HPBW and the hor- 
izontal dashed line the 3 a noise level. 



emission. With these premises, the luminosity is 24 L© 
and the total mass 16 Mq. These values of mass and lu- 
minosity might refer to a possible disk around the YSO 
and should be considered as an upper limit to the true 
disk values since contamination from the dust in the core 
that surrounds the disk might occur. The middle dashed 
curve is a fit to the two sub-mm flux densities of the com- 
pact source (see Sect. I2.2|) . The luminosity is 160 L© and 
the mass 101 Mq. These values should refer to the dens- 
est part of the molecular core. Finally, the upper full line 
is a fit to the two total sub-mm flux densities (compact 
source + halo) and gives a luminosity of 920 L© and a 
mass of 690 Mq. These values are representative of the 
total luminosity of the YSO (or cluster of YSOs) and the 
total mass of the molecular core in which the YSO(s) is 
(are) harboured. 

There are two other ways to derive the mass of the 
molecular core: the first (virial mass) is from the width 
of the line (from gaussian fits) and the source size. We 
used the relation of MacLaren et al. (1988) for a homo- 
geneous sphere. The second uses the intensity of optically 
thin lines and assumes LTE and a given molecular abun- 
dance with respect to H2. The derived masses are given 
in Table |3 and Table respectively. The values obtained 
from HCO + which refer to a resolved source of 0.1 pc in 
size, are similar and consistent, within the approximation 
of the present estimates, with the mass of the dense core 
derived from the fit to the JCMT flux densities. Those 
derived from C 34 S (similar for the two methods) refer to 
a smaller and probably denser central component and are 
consistent with the disk mass obtained from the fit to the 
mm flux densities. Finally, the CH3CN estimates are not 
fully consistent with each other, very likely due to the 
large uncertainty in the molecular abundance assumed. 
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Fig. 14. The SED of the region around the water maser. 
The two black squares are the total 3.3 and 1.2 mm flux 
densities, the filled triangles arc the total flux density and 
the empty triangles the flux density of the compact core 
at 850 and 450 fim, the four circles are the IRAS flux 
densities, the four asterisks the MSX flux densities. The 
three lines denote grey-body fits to the mm and sub-mm 
data discussed in the text. 
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Fig. 15. Overlay of the 850 /im map of the region around 
S235C (contours) with the 21 fim map from MSX (grey 
scale) . The ellipse schematically indicates the position and 
size of the 6 cm radio continuum emission from the S235C 
Hll region (from Israel & Felli I1978j) . 
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Fig. 16. Boltzman plot for the CH 3 CN(5-4), ^=0,1,2,3 
lines. The column densities are source averaged. The 
straight line represents a least squares fit to the data. 

Table 3. Diameter, velocity width and virial masses 



Line 


Diameter 


AV 


Virial Mass 




(pc) 


(km s" 1 ) 


(Mq) 


HCO+(C17) a 


0.1 


2.4 


58 


HCO+(C19) a 


0.08 


3.1 


77 


C 34ga 


0.03 


1.49 


7 


CH 3 CN 


0.06 


2.85 


53 


a Line wings were excluded from fits. 




Table 4. Diameter and mass of the molecular core 


Tracer 


Diameter 


Mass 






(PC) 


(M Q ) 




HCO+ (C17) a 


0.1 


31 




HCO+ (C19) a 


0.08 


10 




C 34 s a 


0.03 


7 




CH 3 CN a 


0.06 


2.3 




mm continuum 


0.054 


16 




a T cx =30 K, 


HCO+]/[H 2 ] 


= lO" 9 , 


[C 34 S]/[H 2 ] = 


[CH 3 CN]/[H 2 ] = 


5 x lO" 10 . 







10" 



3.2. The YSOs 

In this section we examine the possible implications that 
can be derived for the YSOs predicted to exist in the 
molecular core from the NIR observations of Felli et 

ai. mm . 

The NIR source closest to the water maser is Ml (Felli 
et al. H997J) and it is also the reddest one in the entire field. 
Ml is the stellar object with the largest near-IR excess in 
the region and it was proposed by Felli et al. ( 195JZJ) as the 
powering source of the water maser. As shown in Fig.[I]it is 
located ~ 5" SSE of the maser position and only ~ 2" east 
of C19 and it lies projected on the thermal jet. However, 
its distance from the C17 core (~ 10") is greater than the 
possible uncertainty in the near-IR astrometry and makes 
a direct correspondence of Ml with the YSO in the centre 
of C17 very unlikely. Ml could still be associated with 
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the YSO in C17 if it is due to reflected radiation from the 
embedded star, as has been suggested to be the case for the 
(proto)star IRAS 20126+4104 (Cesaroni et al. ITM7|> and 
W75N (Moore et al lTMTl Hunter et al. 11351 Shepherd et 
al. I2003fl . where the NIR emission is not coincident with 
the location of the YSO. This would imply the existence 
of a low density channel with lower extinction through 
which NIR photons emitted by the YSO can escape from 
the dense core. 

We can check the nondetection of a K band source in 
the C17 core by estimating the required stellar extinction 
Ay and consequent column density, and comparing it with 
the column density derived from our mm observations. 
An appropriate lower limit to the K magnitude in the 
region around Ml is 15. To derive Ay we shall assume 
the pre-main sequence (PMS) star to have at least a mass 
of 6 Mq (see Sects. EHfl and 14^1. Using the evolutionary 
tracks of Palla & Stahlcr (JT5J5J, a 5 x 10 5 yr old PMS 
star of such mass at 1.8 kpc would have K=10.9 mag, 
yielding Ay=37 mag. The corresponding column density 
is N(H2) = 3.6 x 10 22 cm~ 2 . The average column density 
for C17 is 3 x 10 22 cm~ 2 from HCO+ and 3.4 x 10 22 cm" 2 
from the mm dust emission. These values are mutually 
consistent and confirm that a 6 Mq YSO hidden inside 
the core could not be observed in the K band. 

Given the proximity of Ml to C19, there is still another 
possibility to consider, namely that Ml is instead coinci- 
dent with a different YSO located in C19 and powering 
the thermal jet. The column density in C19 is lower than 
in C17 (1.4 x 10 22 cm -2 from HCO+), but could even be 
overestimated since no thermal mm continuum is found at 
the same position. This would allow one to detect a lower 
mass star, considering that the putative YSO in C19 must 
have a lower luminosity, but no quantitative statement can 
be made at this stage. In conclusion, it cannot be excluded 
that Ml may be the powering source of the thermal jet and 
the NIR counterpart of a YSO in C19. Higher resolution, 
deeper NIR observations with more accurate astrometry 
are needed to settle this issue. 

3.3. The molecular outflows 

The map of the high velocity HCO + gas (Fig. [SJ shows 
that this is distributed in two outflows, both of them 
originating approximately from the position of the water 
maser. The longest, better denned and more collimated 
one is in the NE-SW direction, with little overlap between 
the red and blue lobes. The collimation is rather high; us- 
ing the water maser as centre position and the HPW of 
the lobes transverse to the maximum length we obtain an 
opening angle of the lobes of ~15°. 

The other outflow is in the NNW-SSE direction, it 
is slightly less extended, less intense and less collimated. 
Its southern blue lobe coincides with the thermal jet seen 
at 3.3 mm and with the H2CS tail, as shown in Fig. Hlsl 
Sulphur-bearing molecules are expected (and found) in 
shocked gas within outflows (Bachiller et al. I2001|) . so the 
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Fig. 17. Overlay of the H2 emission (grey scale from 
Felli et al. IT5571 with the HCO+ outflow (contours). 
The dotted-line contours delineate the 3.3 mm contin- 
uum emission from S235A. The star marks the position 
of S235B, the cross that of the water maser. 

presence of H2CS along the NNW-SSE outflow confirms 
that the molecular tail is part of the outflow. 

Collimated outflows are often associated with shock 
excited H2 emission (see e.g. Cesaroni et al. 1999a). A 
comparison of Fig. [5] with the H2 emission map from Felli 
et al. ( 1 9!]T| is shown in Fig. El Within the uncertain- 
ties of the respective astrometry, the northern blue lobe 
of the NE-SW outflow seems to coincide with a distinct 
feature of the H2 emission that clearly deviates from the 
H2 half ring that surrounds the S235A H11 region. There 
are three other possible associations in Fig.ll7lbetween H2 
features and the bipolar outflows, even though they should 
be checked with more accurate astrometry and more sensi- 
tive H 2 observations: in the northern part, the red lobe of 
the NNW-SSE outflow seems to end up in a bright feature 
of the H2 emission and, in the south, two weaker H2 spots 
are located at the end of the blue lobe of the NNW-SSE 
outflow and of the red lobe of the NE-SW outflow. 

The possible implication is that part of the H2 (the ring 
around S235A) could be excited in the photo dissociation 
region at the interface between the Hn region and the 
molecular cloud, while the H2 associated with the outflows 
could be shock excited. Near-IR spectroscopy of the H 2 
features is needed to settle this point. 

The parameters of the two outflows were derived as- 
suming [HCO + ]/[H 2 ] = 10~ 9 , an excitation temperature 
of 30 K and optically thin emission, and are given in 
Table El The main uncertainty in these values comes 
from the artificial separation of the two outflows in 
the regions where they overlap. All the derived param- 
eters are smaller than those of the outflow found in 
IRAS 20126+4140 (Cesaroni et al. 11557) . as expected con- 
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Fig. 18. Top panel: overlay of the HCO + outflow (con- 
tours) with the 3.3 mm continuum (grey scale). Bottom 
panel: overlay of the HCO + outflow (contours) with the 
H2CS map (grey scale). The cross marks the position of 
the water maser. 



sidering that IRAS 20126+4140 has a bolometric luminos- 
ity about an order of magnitude greater. 

The maximum expansion velocities of the two lobes of 
the NE-SW outflow are not equal. The blue lobe velocity 
reaches up to ~15 km s _1 from the molecular core peak 
velocity, while that of the red lobe is at least twice as 
high. This explains the larger dynamical time scale found 
for the blue lobe, given the fact that the two lobes have 
almost equal sizes. 



3.4. The two velocity peaks 

To study the significance of the two velocity peaks dis- 
played in the HCO + and CH3CN spectra and outlined in 
Sects l2~T~2l and l2~l~3l we have overlaid the HCO+ C17 and 
C19 cores with the 3.3 mm continuum emission (Fig. I19(l 
and the HCO+ outflow (Fig. HI . 



C17 is associated with the 3.3 mm core, the water 
maser, and the centre of the NE-SW outflow. There is 
little doubt that the NE-SW outflow originates from the 
YSO hidden within the 3.3 mm core. 

C19 instead does not have an associated 3.3 mm core 
and is located along the axis of the thermal jet observed 
at 3.3 mm and the blue lobe of the NNW-SSE outflow. At 
the same position and in the same velocity interval (from 
-18 to -20 km s^ 1 ) one sees the jet observed in the H2CS 
line. 

The interpretation of this complex morphology is far 
from clear. The most obvious one is that the molecular 
core is composed of two distinct components with different 
velocities. However, in this case there would be a major 
difference between C17 and C19 due to the lack of mm 
continuum emission from C19: either the temperature of 
the dust in C19 is much lower than in C17 or the abun- 
dance of HCO + and CH3CN in C19 is much higher than in 
C17. If this were the case it would be of interest to know if 
the thermal jet and the NNW-SSE outflow originate from 
C19, but this cannot be derived from the present data. 

An alternative explanation might be that the accelera- 
tion of the molecular gas from its interaction between the 
thermal jet and the molecular outflow is capable to alter in 
this area the velocity observed in the HCO + and CH3CN 
lines, which would be the weighted mean of the gas en- 
trained and the quiescent surrounding molecular gas. In 
IRAS 18162-2048 Benedettini et al. (2003) have found 
that the momentum deposited by the outflow can break 
the molecular cloud apart, shifting part of the emission 
towards bluer or redder velocities. If all the momentum in 
the blue lobe of the NNW-SSE outflow (13 M©km s" 1 ) 
is deposited in the molecular cloud, given the velocity off- 
set between C17 and C19 (2 km s _1 ), up to 6 M can be 
moved. This mass is not too different from that of C19 
(see Table 0} and, if accelerated, would definitely create a 
secondary velocity peak in the molecular line emission. 

3.5. The thermal jet 

The detection of thermal jets in YSOs is by no means a 
new discovery since many have been found in low lumi- 
nosity YSOs (Anglada et al. I1998|l . Usually they are as- 
sociated with molecular outflows, even though on a much 
smaller scale, and are elongated along the axis of the out- 
flow. The origin of the ionization of the thermal jet is 
not photoionization, because at the low bolometric lu- 
minosities the implied stellar temperatures are too low. 
The alternative mechanism proposed is ionization from 
UV photons produced in the shock of the neutral wind 
from the YSO against the surrounding high density gas 
(Curiel etal. H2H3 H2H2I Torrelles et a. H555|l . Anglada 
et al. I|1998|l have found an empirical relation between 
the momentum rate in the molecular outflow (P) and 
the radio continuum luminosity of the jet at centimeter 
wavelengths (S v d 2 ). Using their Eq. (3) and the extrapo- 
lated centimeter flux density of the thermal jet, we derive 
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Fig. 19. Overlay of the HCO+ C17 (top) and C19 (bot- 
tom) cores (contours) with the 3.3 mm map (grey scale). 
The cross marks the position of the water maser. 
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Fig. 20. Overlay of the HCO + outflow (contours) with 
C17 (top panel) and C19 (bottom panel) HCO + cores 
(grey scale). The cross marks the position of the water 
maser. 



P = 9 10 3 Mq km s x . This value can be compared to 
that derived for the blue lobe of the NNW-SSE outflow 
using the values of Table i.e. P = 10~ 3 M km s _1 . 
Although the two values differ by almost a factor 10, we 
believe that they are in reasonable agreement within the 
uncertainties. In fact, on the one hand the relation by 
Anglada et al. fUSEI has a spread of about an order of 
magnitude, on the other the masses (and hence the mo- 
menta) in Table[3]depend on the HCO + abundance, which 
may differ significantly from the value adopted by us. 



3.6. The disk 

Highly collimated bipolar outflows have been found to 
be associated with disks in several YSOs, as for instance 
in HH 30 (Bur rows e t al. ITMfl and IRAS 201026+4140 
(Cesaroni et al 119971 ITMIafl . The detected disks rotate 
around the central (proto)star and lie in a plane perpen- 
dicular to the bipolar outflows. One way to detect the disk 
is by revealing the change of the LSR velocity along a di- 
rection perpendicular to the bipolar outflow. To search for 
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the possible presence of a disk in the centre of the HCO + 
outflow we have used the highest resolution available in 
our data and the strongest molecular line, i.e. the C 34 S 
observations. 

As pointed out in Sect. 12.1.41 from the comparison 
of single dish and interferometric observations, the C 34 S 
emission is dominated by the quiescent gas of the ex- 
tended molecular cloud. Therefore, in order to avoid con- 
fusion we have selected the extreme velocity ranges where 
the emission is still above the noise level, thus averaging 
the emission from -20.86 to -18.9 km s _1 (blue side) and 
from -15.81 to -13.86 km s _1 (red side). The flux density 
contained in these intervals is ~17% of that within the 
above velocity limits. The resulting maps are shown in 
Fig- El Two unresolved components are found (Fig. 121b,) 
very close to the water maser and to the 1.2 mm core 
(Fig. 121b - ). which also coincide with the centre of the NE- 
SW outflow (Fig. 121b). The red and blue peaks are off- 
set by ~3.5" and the line that connects them is almost 
perpendicular to that of the NE-SW outflow. This might 
suggest that one is observing a disk rotating around the 
YSO and perpendicular to the outflow axis. One could 
possibly object that the two unresolved components rep- 
resent the feet of the NNW-SSE outflow, given also the 
good correspondence in velocity and position. While we 
cannot exclude this hypothesis, we do not feel this is the 
case because the C 34 S molecule is a good tracer of high- 
density gas. To confirm our hypothesis we have made a 
velocity-position plot along the direction of the plane of 
the disk. This is shown in Fig. which clearly indicates 
the existence of a velocity gradient. 

Assuming that we are dealing with a rotating disk, 
it is possible to estimate the mass of the star+disk sys- 
tem from the measured radius and corresponding rota- 
tion velocity and check if this value is consistent with the 
other indications. The radius can be computed from the 
separation between the red- and blue-shifted C 34 S gas in 
Fig. I21h. which at the assumed distance of 1.8 kpc corre- 
sponds to 3150 AU. The rotation velocity is taken equal 
to the maximum velocity offset between the most blue- 
and red-shifted channels in which the C 34 S emission is 
still above the noise, and turns out to be ~5 km s _1 . 
The inferred dynamical mass of the disk+star system is 
>22 Mq. Although this is not a demonstration of the ex- 
istence of a disk, such a value compares well to the disk 
mass found in other YSOs with similar mass, as for in- 
stance IRAS 20126+4104 (Cesaroni et al.ESZ). 

The lower limit comes from the unknown inclination 
angle of the disk with respect to the line of sight. This 
value is larger than the disk mass derived from the fit to 
the mm flux densities. Given all the uncertainties, this 
result is consistent with the hypothesis of a disk of mass 
~16 Mq, as derived from the fit to the SED (see Sect. 13. l|) . 
rotating around a YSO of >6 Mq. 



4. The nature of the YSOs in the molecular core 

The outflows observed in the HCO + line witness the exis- 
tence of embedded YSOs powering them. The question 
we want to address here is whether these are low- or 
high-mass stars. The molecular core is significantly colder 
(30 K) than a "hot core" , whose temperature is typically 
above 100 K, and also the line width is smaller than that 
found in hot cores: since hot cores are believed to host 
newly born massive stars (see Kurtz et al. I2000|) . such a 
temperature and line width difference might suggest that 
one is observing lower mass YSOs in the molecular core. 
To find out whether this is the case, one has to obtain a 
better constrained estimate of the luminosity of the YSOs. 

The presence of two outflows guarantees that at least 
two distinct YSOs must be present in the core. The possi- 
bility that multiple outflows may originate from the same 
YSO, although it cannot be excluded a priori, seems less 
supported by the orthogonal position of the two outflows. 

According to Churchwell (1997), the luminosity of the 
star powering the outflow can be inferred from the mass 
loss rate. Using his Eq. (1) and our results from Table 03 
one obtains 3400 Lq and 170 Lq for the YSOs powering 
respectively the NE-SW and NNW-SSE outflows. Clearly, 
the second cannot contribute significantly to the total lu- 
minosity estimated in Sect. l3.il which is dominated by the 
most massive YSO. 

Another luminosity estimate can be calculated from 
the temperature obtained from the CH3CN lines. We mea- 
sure 30 K at a distance of 0.03 pc (the radius of the 
CH3CN emitting region; see TableOJ from the core centre. 
If the gas is heated by a central star, then it is possible 
to compute the stellar luminosity. The result depends on 
the geometry of the core. If this is spherically symmetric, 
one can use the relation of Scoville & Kwan l|1976[) . which 
gives a luminosity of 220 Lq for k v oc v 5 . Much higher 
luminosities are predicted, up to 2700 Lq, if instead we 
use the more canonical value of k v oc v and Eq. (31) of 
Goldreich & Kwan (1974). For a disk-like geometry, in- 
stead, even an 03 star cannot justify a temperature as 
large as 30 K at a distance of 0.03 pc from the star. It 
must be noted that this result does not rule out the possi- 
bility that one is indeed observing a circumstellar disk, as 
discussed in Sect. 13.61 but this must be "active", namely 
heated internally by viscosity. 

Finally, one can infer the luminosity from the mass of 
the star estimated in Sect. 13.61 The latter is >6 Mq, which 
implies >1000 Lq for a zero-age main-sequence star or 
alternatively >3100 Lq for a massive protostar (Behrend 
& Maeder EIII1)- 

All previous luminosity estimates must be compared to 
the value obtained from the SED in Fig.^J which ranges 
from 160 to 900 Lq. However, the uncertainty in this value 
is not negligible, as the SED is known only in the millime- 
ter and sub-millimctcr regions of the spectrum. On the 
other hand, the other luminosity estimates are also af- 
fected by significant errors. In fact, the relationship used 
to derive the luminosity from the mass loss rate has a 
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spread of a factor ~3, while the value obtained from the 
core temperature is proportional to T 5 , which makes it 
very sensitive to small errors on T . In conclusion, we be- 
lieve that the estimates obtained are consistent with each 
other within the uncertainties and that the YSO powering 
the NE-SW flow is likely to be a B (proto)star with lumi- 
nosity in the range 1000-3500 L , while the NNW-SSE 
outflow originates from a less massive (proto)star. 

We note that for this range of luminosities, the mass 
infall rate derived from the relation of Churchwell l|1997f) 
and Behrend & Maeder (|5001 ) turns out to be of the or- 
der of 10~ 3 M & yr" 1 , not too different from the mass infall 
rate derived from the infall velocity. Even though the close 
correspondence between the two numbers should be taken 
with caution, given the large uncertainties in the derived 
quantities and the approximations used, it suggests that 
gravitational energy might be the main source of luminos- 
ity for the YSO powering the water maser. 

5. Conclusions 

The detection of a mm and sub-mm core at the position 
of the water maser in S235A-B, as well as of a core in all 
the molecular lines finally settles the question which is the 
source that provides the energy for the maser emission and 
indicates that a new YSO (or a cluster of YSOs) is located 
in this region and is the powering source of the maser 
emission. The two close-by nebulosities S235A-B do not 
seem to be directly related to the YSO, even though be- 
longing to the same star forming complex. Consequently, 
star formation has not been occurring simultaneously in 
the extended CO molecular cloud, but rather progressing 
slowly from the outer parts of the molecular cloud, where 
S235A and B are placed, to the denser inner part of it, 
where new YSOs are just being formed. 

The present example and a look at the literature on 
this star forming region indicated the complexity in find- 
ing new intermediate luminosity YSOs when they are lo- 
cated in regions where high-mass star formation is already 
in progress and where their emission is completely obliter- 
ated by that of the more evolved regions, especially when 
the search occurs in the classical domains which were orig- 
inally used to find them, e.g. cm radio continuum, near 
IR and Ha emission. Only water masers, as well as other 
masers, e.g. methanol (Walsh et al. 2003J), have remained 
as good indicators of YSOs (when they are present) and 
only high resolution mm continuum and molecular line 
observations are able to pinpoint the YSO emission from 
that of other close-by more evolved sources (Jenness et 
aLHSEJ. 

The detection of two bipolar outflows in S235 is an- 
other good indicator of recent star formation. It indicates 
that more than one YSO might be present or, but more 
difficult to explain with the present data, multiple episodes 
of bipolar outflow might originate from the same YSO. 
Larger scale bipolar outflows detected in this region with 
lower resolution might be remnants of previous episodes 
of star formation, but do not seem to be directly related 



to the two outflows produced by the presently detected 
YSOs. 

Sulphur-bearing molecules (H2CS and SO2) have been 
detected, similarly to hot cores. Part of the emission ap- 
pears related with the outflows, suggesting a shock origin, 
but most of it clearly arises from the molecular core. 

The luminosity of the YSO powering the NE-SW flow 
is estimated to be in the range 1000-3500 Lq, less than 
that of S235A (1.1 10 4 L©) and the dynamical time of the 
order of 10 4 yr, much younger than that of the H11 region 
around S235A (the canonical value for well developed H11 
region is 10 6 yr). The temperature of the molecular gas of 
the core (~30 K) is smaller than that found in hot cores 
(>100 K) further confirming that the newly detected YSO 
is intermediate between the massive objects found in hot 
cores (Kurtz et al. l2000|l and the low mass class objects 
(Andre et al. 1155511 . 

In conclusion, the newly found YSO in between S235A 
and S235B is a rare example of an intermediate luminos- 
ity object and its morphology may be used to connect 
early evolutionary phases of massive stars with those of 
low mass protostars of class 0-1. 

We find weak evidence for a possible disk in the C 34 S 
line, perpendicular to the direction of the main bipolar 
outflow. If confirmed, this is a further proof of how the 
complex morphology of star forming regions in the earliest 
phases can be unveiled only by high resolution molecular 
observations in high density tracers that can separate the 
weak emission of the disk from the overwhelming emission 
of the molecular cloud. 

While revealing a new wealth of information, at the 
same time the latest observations leave some unanswered 
questions, as for instance the origin of the thermal jet and 
the nature of the two-velocity structure of the molecular 
cloud. Are there multiple YSOs in the region of the mm 
core or a single YSO with multiple episodes of outflows? 
Why there is no thermal jet associated with the main bipo- 
lar outflow? Why is the velocity of the water maser fea- 
tures so different from that of the molecular core? Work 
is in progress to try to give an answer to these questions. 
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Fig. 21. a. Blue and red shifted emission of the C 34 S line 
(contour levels from 60 to 180 by 60 mJy/beam). b. Same 
as a overlaid with the 1.2 mm map (grey scale), c. Same 
as a overlaid with the most intense part of the HCO + 
outflow (contours). The cross marks the position of the 
water maser. 
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Fig. 22. Position-velocity plot of the C S along the di- 
rection of the plane of the disk, averaged over a width of 
4" in the perpendicular direction. The dotted lines define 
the velocity interval which has been used to produce the 
maps of Fig. [21 The zero offset position corresponds to 
the water maser. Contours are from 0.05 to 0.55 in steps 
of 0.05 Jy/beam. 



